Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  618 1567 (N24E40 L=168)  2898    551        0 ( 464/1448- 834/1729) -13.8 E69-E20 N19-N32
      2  990 1430 (N22E10 L=198)   267     83        0 ( 957/1411-1034/1453) -11.3 E12-E07 N21-N23
      3 1497 1342 (N27W25 L=233)   384    133        0 (1456/1293-1539/1389) -35.5 W21-W28 N24-N29
      4  863 1239 (N08E14 L=195)   541    211        0 ( 796/1215- 961/1291)   0.0 E18-E08 N06-N10
      5  682 1185 (N02E24 L=185)   419     98        0 ( 636/1164- 736/1201)  18.4 E26-E20 N00-N03
      6 1136 1100 (N05W05 L=214)   438    205        0 (1112/1007-1156/1196) -79.7 W03-W06 N00-N11
      7  559  826 (S22E26 L=182)  1061    294        1 ( 506/ 733- 657/ 931) -42.3 E30-E18 S26-S15
      8 1445  371 (S33W50 L=258)   695    103        2 (1396/ 352-1496/ 395)  14.0 W45-W53 S34-S30
  

    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