Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  805 1339 (N12E20 L=103)   319     88        0 ( 774/1313- 851/1373)  -9.5 E22-E16 N10-N13
      2  692 1141 (S02E22 L=101)   486    161        0 ( 641/1069- 750/1217) -29.1 E26-E17 S04-N01
      3  940 1078 (S00E06 L=117)   547    221        1 ( 911/ 981- 968/1184) -80.5 E08-E04 S06-N06
      4  266 1099 (S12E53 L= 70)   435     87        1 ( 250/1059- 277/1148) -59.0 E55-E51 S14-S09
      5  455 1053 (S11E36 L= 87)   284     74        1 ( 424/1037- 490/1076)  36.9 E38-E34 S12-S09
      6  668  877 (S18E19 L=104)  1051    413        1 ( 488/ 772- 788/1012)  51.3 E30-E13 S27-S07
      7 1866  985 (N14W61 L=185)   469    200        0 (1826/ 971-1913/ 999)   9.5 W56-W68 N12-N15
      8 1409  877 (S02W26 L=149)   256     63        0 (1394/ 843-1424/ 912) -63.4 W24-W27 S04-N00
      9  645  739 (S27E18 L=105)   236     98        0 ( 612/ 695- 681/ 774)  60.3 E20-E16 S30-S23
     10 1335  580 (S21W30 L=153)   261     58        0 (1303/ 554-1358/ 610)  63.4 W28-W31 S23-S19
     11 1718  531 (S16W63 L=186)   537     60        0 (1690/ 503-1745/ 554)  45.0 W61-W65 S17-S13
     12 1216  363 (S38W32 L=155)  1092    243        1 (1092/ 342-1333/ 382)  23.6 W22-W39 S41-S34
  

    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