Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1263 1475 (N33W14 L= 76)   351     62        1 (1216/1455-1295/1490)   9.5 W10-W16 N31-N33
      2 1410 1407 (N29W24 L= 87)   564    108        2 (1375/1372-1455/1457) -51.3 W21-W27 N26-N32
      3  421 1281 (N13E43 L= 19)   875    184        2 ( 394/1216- 456/1342) -63.4 E46-E40 N08-N16
      4  225 1220 (N07E62 L=  1)   957    168        2 ( 197/1158- 252/1270)  74.1 E65-E58 N03-N10
      5 1192  831 (S09W12 L= 74)   432    144        1 (1171/ 759-1214/ 898) -77.5 W10-W13 S13-S04
      6  911  740 (S17E05 L= 57)   265    115        0 ( 877/ 695- 952/ 782)  56.3 E07-E02 S19-S13
      7 1613  572 (S23W47 L=109)   459    106        0 (1597/ 527-1634/ 610) -51.3 W44-W48 S25-S20
      8 1532  451 (S32W44 L=107)   461     97        1 (1491/ 431-1573/ 481)  15.9 W40-W47 S33-S29
  

    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