Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  827 1731 (N35E36 L=259)   522    120        1 ( 784/1704- 875/1760) -11.3 E40-E30 N33-N36
      2  619 1644 (N24E45 L=250)  1060    265        2 ( 580/1542- 643/1719) -42.5 E51-E38 N17-N28
      3  271 1487 (N07E69 L=225)   611     95        1 ( 256/1454- 289/1513) -26.6 E73-E66 N05-N08
      4  550 1405 (N08E39 L=256)   440     81        2 ( 516/1377- 580/1439)  11.3 E41-E36 N05-N10
      5  526 1259 (S01E36 L=259)   230     88        0 ( 513/1214- 549/1300) -45.0 E37-E33 S03-N01
      6  236 1157 (S13E58 L=237)   748     98        3 ( 216/1119- 258/1198) -38.7 E60-E54 S15-S10
      7  181 1112 (S17E63 L=232)   891     94        2 ( 159/1061- 204/1153)  80.5 E64-E60 S19-S13
      8  659  993 (S14E21 L=274)   406     80        1 ( 626/ 956- 699/1029)  90.0 E22-E18 S16-S11
      9  194  836 (S33E57 L=237)   497     76        1 ( 172/ 815- 219/ 858) -11.3 E59-E54 S34-S32
     10  628  754 (S28E18 L=276)   314     60        1 ( 602/ 725- 656/ 777)  53.1 E19-E16 S30-S26
     11  478  619 (S42E29 L=266)   635    111        1 ( 437/ 580- 517/ 659)  45.0 E32-E25 S45-S38
     12 1560  365 (S27W59 L=354)   451    100        0 (1528/ 351-1590/ 384)   8.1 W55-W62 S27-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