Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  841 1838 (N55E20 L=162)   609     84        1 ( 794/1810- 884/1856)   5.7 E25-E15 N51-N57
      2  292 1463 (N25E58 L=124)   939    125        3 ( 266/1410- 318/1528) -45.0 E61-E52 N21-N29
      3 1626 1508 (N28W45 L=227)   409     82        0 (1585/1479-1660/1528) -26.6 W41-W47 N26-N29
      4  451 1420 (N22E40 L=142)   549    145        2 ( 412/1368- 485/1482) -41.6 E45-E36 N18-N26
      5  925 1322 (N15E06 L=176)   222     56        0 ( 894/1310- 964/1334) -14.0 E08-E03 N13-N15
      6  378 1240 (N10E43 L=139)   411     81        1 ( 348/1203- 398/1276)  63.4 E45-E41 N08-N12
      7  495 1210 (N08E34 L=148)   228     44        0 ( 474/1192- 519/1229)  33.7 E35-E32 N07-N09
      8  451 1115 (N03E37 L=145)  1254    294        2 ( 400/ 985- 534/1225)  77.0 E40-E30 S04-N09
      9 1109  833 (S15W05 L=187)  1057    322        3 (1025/ 703-1158/ 975)  66.4 W00-W08 S22-S06
     10 1625  785 (S17W41 L=223)   397    119        0 (1572/ 767-1671/ 810)   7.1 W36-W44 S17-S15
     11 1096  474 (S38W05 L=187)  1226    393        2 ( 919/ 445-1291/ 500)   6.1 E08-W20 S40-S36
     12  983  304 (S52E04 L=178)   536    121        0 ( 921/ 288-1041/ 322)   0.0 E09-W01 S53-S50
  

    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