Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1578 1372 (N12W44 L=204)   389     87        0 (1549/1333-1605/1413) -90.0 W42-W45 N09-N15
      2 1889 1174 (S08W67 L=226)  1159    135        2 (1864/1134-1912/1211)  78.7 W63-W70 S10-S05
      3 1774  967 (S16W48 L=207)   457    122        1 (1743/ 909-1801/1025)  56.3 W44-W51 S19-S13
      4 1672  985 (S12W39 L=199)   350    124        2 (1641/ 948-1701/1024)  76.0 W37-W41 S14-S09
      5 1374  921 (S07W18 L=177)   823    284        4 (1281/ 802-1438/1006) -75.1 W13-W20 S15-N00
      6  162  886 (N14E67 L= 92)  1414    206        6 ( 125/ 839- 198/ 941) -50.2 E71-E62 N11-N17
      7 1259  507 (S30W02 L=161)  1170    441        4 (1108/ 381-1420/ 662)  22.9 E11-W15 S35-S23
      8  436  469 (S14E55 L=105)  2996    539        9 ( 365/ 387- 555/ 557) -58.4 E62-E48 S21-S07
      9  354  443 (S15E65 L= 94)   771    185        4 ( 322/ 418- 382/ 466) -76.0 E71-E62 S17-S13
  

    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