Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1683 (N40E43 L=123)   651    120        1 ( 595/1635- 674/1725) -24.8 E52-E39 N36-N42
      2  424 1613 (N32E61 L=106)   618     75        0 ( 405/1573- 445/1654) -24.4 E66-E55 N29-N34
      3 1265 1588 (N41W12 L=179)   770    157        1 (1203/1553-1338/1632) -14.0 W06-W18 N38-N44
      4  637 1435 (N23E32 L=135)   654    140        2 ( 601/1378- 673/1501)  90.0 E33-E29 N18-N27
      5 1487 1224 (N18W29 L=195)   333    127        0 (1464/1172-1526/1279) -56.3 W26-W31 N15-N21
      6  285  887 (S15E52 L=115)  4018    509        3 ( 186/ 661- 386/1052)  40.3 E76-E42 S31-S03
      7 1391  912 (S02W24 L=191)   377    183        0 (1312/ 873-1451/ 930)  -5.7 W18-W28 S03-N00
      8 1688  730 (S10W50 L=217)   445     80        0 (1676/ 684-1700/ 771) -71.6 W48-W51 S13-S07
      9 1224  748 (S14W15 L=182)   329     95        1 (1180/ 733-1273/ 763)   0.0 W12-W18 S14-S12
     10 1210  558 (S26W18 L=185)   260     90        0 (1168/ 537-1252/ 575)  31.0 W15-W20 S28-S24
     11 1632  508 (S26W55 L=221)  1233    155        1 (1563/ 474-1694/ 542)  24.4 W49-W60 S28-S23
     12 1249  438 (S34W24 L=191)   670    136        1 (1183/ 420-1304/ 458)   0.0 W18-W28 S36-S32
  

    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