Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1068 1811 (N52W13 L= 66)   557    150        0 (1025/1766-1115/1853)  30.3 W08-W20 N48-N56
      2 1193 1669 (N38W19 L= 71)   240     48        0 (1158/1655-1220/1686)   9.5 W15-W21 N37-N39
      3 1627 1473 (N22W46 L= 98)   279     54        0 (1597/1458-1654/1491) -26.6 W43-W47 N20-N22
      4  585 1340 (N21E26 L= 26)   710    221        1 ( 517/1285- 672/1398) -39.3 E31-E20 N16-N25
      5 1405 1054 (S03W23 L= 75)   193     87        0 (1384/1008-1424/1104)  68.2 W21-W24 S05-N00
      6  152 1005 (N04E65 L=348)  1714    217        4 ( 118/ 895- 186/1069) -68.2 E69-E60 S02-N08
      7  226 1025 (N05E56 L=357)   340     85        1 ( 209/ 982- 246/1068) -78.7 E57-E53 N02-N07
      8 1600 1010 (S06W36 L= 89)   237     34        1 (1582/ 981-1621/1033) -90.0 W35-W38 S07-S04
      9  669  341 (S43E37 L= 16)   664    126        1 ( 620/ 324- 737/ 363)   5.2 E40-E29 S44-S40
  

    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