Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1577 (N33W28 L=316)  1115    269        0 (1063/1550-1291/1617) -13.2 W20-W37 N29-N36
      2 1173 1311 (N18W17 L=304)   194     58        0 (1155/1287-1200/1331)  14.0 W14-W18 N17-N19
      3  545 1305 (N35E24 L=264)   253     83        0 ( 511/1272- 582/1326)   9.5 E26-E20 N33-N35
      4  928 1244 (N21W00 L=288)   356     67        0 ( 901/1209- 949/1285)  63.4 E01-W02 N19-N23
      5  722 1241 (N27E12 L=275)   254     34        0 ( 696/1224- 751/1257)   0.0 E14-E09 N25-N27
      6  323  639 (N01E58 L=230)   399     69        0 ( 284/ 627- 367/ 649)  -6.3 E62-E53 N00-N01
      7  345  522 (S07E62 L=225)   928    116        5 ( 290/ 492- 382/ 550) -24.0 E67-E58 S08-S04
      8  480  371 (S19E60 L=228)   497     45        0 ( 463/ 339- 498/ 395)  68.2 E62-E56 S21-S16
      9  607  321 (S25E52 L=235)   419     56        1 ( 586/ 295- 629/ 353)  45.0 E54-E49 S26-S22
  

    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