Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  778 1387 (N33E05 L= 96)   351     79        1 ( 737/1363- 823/1406)  -9.5 E07-E01 N32-N34
      2 1857 1198 (S10W61 L=162)  1023    213        5 (1825/1164-1899/1245) -56.3 W56-W66 S12-S06
      3 1530 1108 (S03W31 L=132)   959    230        3 (1495/1005-1575/1232)  72.9 W28-W35 S10-N03
      4  872 1192 (N20E04 L= 96)   281     99        1 ( 850/1140- 901/1229)  38.7 E06-E01 N16-N21
      5 1813 1140 (S11W54 L=155)   387     73        1 (1792/1115-1840/1161)   0.0 W51-W57 S12-S09
      6 1861 1043 (S18W56 L=157)   577    119        2 (1833/1013-1893/1080)  15.9 W52-W60 S19-S16
      7  304  943 (N19E49 L= 52)   464    149        0 ( 258/ 904- 360/ 988) -45.0 E52-E44 N15-N22
      8 1750  803 (S28W41 L=141)  1028    272        1 (1717/ 676-1790/ 923)  63.4 W38-W46 S36-S21
      9 1220  655 (S20W01 L=101)  1072    388        2 (1163/ 511-1265/ 869)  55.6 E06-W08 S27-S08
     10  331  611 (S01E58 L= 42)   647    122        1 ( 316/ 561- 350/ 660)  33.7 E61-E54 S03-N01
     11  470  490 (S11E51 L= 49)   342    163        0 ( 415/ 469- 519/ 517)  -5.7 E56-E46 S11-S09
     12  970  390 (S29E23 L= 77)  1009    234        0 ( 874/ 348-1037/ 467)  -4.8 E28-E15 S32-S25
  

    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