Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  799 1676 (N51W06 L=264)   333     70        0 ( 781/1641- 816/1711)  56.3 W04-W08 N47-N53
      2 1464 1687 (N29W54 L=312)   351     42        0 (1445/1664-1485/1710)   8.1 W50-W57 N28-N30
      3  407 1261 (N33E36 L=222)   263    109        1 ( 358/1238- 459/1289) -24.0 E40-E31 N31-N36
      4  310 1195 (N31E46 L=213)   679    119        1 ( 257/1158- 361/1229)  10.3 E51-E40 N29-N32
      5  449 1096 (N22E34 L=224)   319     78        0 ( 407/1080- 496/1115) -29.7 E37-E30 N20-N24
      6 1703  712 (S33W38 L=296)   402     53        1 (1676/ 688-1736/ 738)   0.0 W35-W41 S33-S31
      7 1418  653 (S27W14 L=272)   283     79        1 (1378/ 624-1445/ 685)   0.0 W10-W16 S28-S25
      8 1709  648 (S37W38 L=297)   343     42        0 (1688/ 623-1731/ 670)  26.6 W36-W41 S38-S35
      9  529  423 (S18E51 L=208)   330     95        1 ( 504/ 381- 551/ 449)  20.6 E55-E47 S20-S17
     10 1457  293 (S58W11 L=269)   407     42        0 (1429/ 278-1480/ 313) -15.9 W07-W14 S60-S55
     11 1231  184 (S59E18 L=240)  1546    176        3 (1164/ 145-1319/ 234) -17.4 E27-E08 S62-S53
     12 1097  155 (S56E33 L=225)   567     56        1 (1072/ 137-1121/ 177)  78.7 E37-E29 S58-S53
  

    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