Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1188 (N00E09 L=188)   272    161        0 ( 902/1138- 992/1243)  68.2 E11-E06 S01-N04
      2  680 1035 (S15E20 L=178)   615    106        2 ( 612/ 992- 735/1094)  58.0 E23-E17 S18-S09
      3 1820  843 (N09W58 L=256)  1751    281        6 (1782/ 754-1865/ 926) -33.7 W52-W68 N05-N13
      4  800  845 (S22E08 L=190)   275     95        0 ( 767/ 822- 849/ 878)   0.0 E10-E04 S23-S20
      5 1393  777 (S09W28 L=226)   752    150        4 (1372/ 694-1410/ 847) -60.9 W24-W31 S13-S04
      6 1001  764 (S21W06 L=204)   260     96        1 ( 950/ 750-1053/ 778)  18.4 W03-W09 S22-S19
      7 1349  707 (S14W28 L=226)   218     59        1 (1329/ 676-1368/ 730) -33.7 W26-W29 S15-S12
      8  576  649 (S40E19 L=179)   826    238        3 ( 478/ 631- 673/ 678)  17.4 E27-E11 S43-S37
      9 1151  501 (S31W25 L=223)   393    126        2 (1092/ 465-1208/ 536)   0.0 W20-W30 S32-S30
     10 1498  459 (S22W49 L=247)  2231    739        4 (1262/ 409-1713/ 509)  18.9 W33-W74 S30-S15
     11 1216  261 (S43W46 L=243)   425     73        0 (1185/ 236-1240/ 281)   0.0 W42-W49 S44-S42
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y