Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1118 1874 (N58E30 L=219)  1007     84        1 (1094/1832-1144/1918) -42.9 E39-E24 N52-N65
      2  762 1735 (N34E42 L=207)   404     69        0 ( 747/1699- 780/1772) -35.5 E45-E38 N31-N36
      3 1291 1653 (N43E00 L=249)   556     71        0 (1265/1622-1317/1688) -31.0 E03-W02 N40-N45
      4 1593 1584 (N49W30 L=278)   865     91        2 (1550/1553-1636/1607)   0.0 W24-W34 N46-N51
      5 1685 1495 (N45W40 L=289)   416     55        1 (1669/1467-1709/1527) -90.0 W37-W42 N42-N47
      6 1637 1250 (N26W34 L=283)  2114    503        4 (1492/1056-1774/1415) -25.8 W20-W49 N17-N32
      7  184  996 (S24E63 L=186)  1568    160        3 ( 148/ 906- 224/1067)  48.0 E69-E57 S29-S19
      8 1079  592 (S27W15 L=264)  1654    393        4 ( 912/ 513-1249/ 646)  45.0 W06-W24 S36-S19
      9 1330  524 (S24W34 L=283)   415    160        0 (1273/ 476-1377/ 584)  77.5 W32-W35 S28-S19
     10 1433  342 (S32W54 L=303)  2287    352        2 (1258/ 319-1561/ 374)  22.6 W39-W65 S38-S26
     11  900  227 (S60W25 L=274)   372     52        1 ( 874/ 204- 927/ 250)  63.4 W21-W28 S62-S56
     12  980  160 (S63W44 L=293)  1712    158        3 ( 946/ 102-1029/ 221)  -8.1 W31-W73 S66-S56
  

    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