Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1601 (N51E12 L=264)   407     52        0 ( 618/1568- 646/1638)  51.3 E13-E09 N48-N54
      2 1297 1578 (N28W35 L=311)   420     60        0 (1260/1560-1338/1594) -18.4 W32-W38 N27-N29
      3  392 1159 (N28E38 L=237)   484     86        0 ( 355/1119- 423/1198) -50.2 E41-E35 N24-N30
      4  475 1074 (N21E33 L=243)   354     89        0 ( 428/1057- 513/1093) -15.9 E36-E29 N19-N21
      5  196  928 (N18E62 L=214)   505     53        1 ( 170/ 900- 218/ 955)  45.0 E65-E58 N16-N19
      6 1849  791 (S33W56 L=332)  1655    119        1 (1814/ 729-1880/ 846)  45.0 W51-W61 S36-S29
      7 1180  544 (S27E04 L=272)   345    122        0 (1129/ 498-1226/ 582)  14.0 E08-W00 S28-S25
      8 1325  221 (S60E07 L=269)   802     76        2 (1289/ 196-1369/ 245) -45.0 E11-E01 S63-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