Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1648 (N27E33 L=147)   329     73        1 ( 712/1602- 746/1681) -56.3 E34-E30 N23-N29
      2  345 1629 (N21E68 L=111)   652     92        2 ( 317/1592- 380/1649) -14.9 E75-E60 N18-N22
      3  373 1531 (N14E56 L=123)  2048    226        9 ( 317/1435- 429/1594)  79.7 E65-E50 N07-N18
      4 1780 1284 (N25W48 L=227)  1905    382        1 (1662/1150-1926/1378) -10.6 W36-W68 N21-N28
      5  302 1122 (S11E50 L=129)   348     78        0 ( 275/1098- 335/1152)  31.0 E52-E47 S12-S09
      6  262 1023 (S17E53 L=126)  1164    199        3 ( 212/ 944- 306/1099) -71.6 E58-E48 S22-S13
      7  596  751 (S30E24 L=156)   254     64        1 ( 569/ 730- 616/ 780)  45.0 E25-E21 S31-S27
      8  725  737 (S29E14 L=165)   604    183        1 ( 638/ 711- 817/ 772)   0.0 E20-E07 S29-S27
      9 1827  597 (S13W73 L=252)   744    109        2 (1808/ 542-1844/ 643) -56.3 W68-W76 S15-S09
     10  279  599 (S44E59 L=120)   730     79        2 ( 235/ 574- 320/ 628)  10.0 E69-E52 S46-S41
     11  431  424 (S56E48 L=131)  1009    125        1 ( 377/ 386- 488/ 463)  12.5 E64-E37 S59-S52
     12 1615  383 (S31W66 L=245)  3206    272        5 (1482/ 339-1739/ 439)  27.6 W53-W80 S36-S24
  

    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