Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1090 1856 (N51E23 L=195)   382     46        0 (1063/1833-1115/1876) -15.9 E26-E19 N48-N52
      2 1258 1353 (N17W06 L=225)  1098    371        1 (1189/1266-1341/1480) -39.8 W00-W12 N13-N23
      3  532 1431 (N07E40 L=179)   737    149        2 ( 478/1381- 590/1476)  49.4 E42-E36 N02-N10
      4  691 1418 (N09E29 L=190)   504    162        0 ( 628/1358- 750/1454) -11.3 E33-E23 N06-N09
      5 1667 1167 (N17W37 L=256)   954    222        2 (1603/1112-1751/1253) -12.1 W30-W44 N14-N20
      6  241  976 (S24E55 L=164)   416     69        0 ( 228/ 935- 252/1016)  78.7 E56-E53 S26-S21
      7  932  818 (S21E01 L=218)   426    192        0 ( 849/ 797-1025/ 844)  33.7 E05-W04 S23-S17
      8  266  749 (S38E54 L=164)   736    122        1 ( 245/ 668- 282/ 795)  53.1 E60-E51 S43-S35
      9  448  670 (S40E34 L=184)   952    142        1 ( 393/ 621- 489/ 736)  39.3 E41-E30 S44-S35
     10 1736  557 (S14W64 L=283)   862    102        1 (1705/ 510-1761/ 591)  90.0 W59-W68 S17-S11
     11 1557  303 (S33W71 L=290)   688     72        0 (1529/ 273-1583/ 338)   0.0 W63-W81 S33-S31
  

    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