Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1275 1475 (N33W15 L= 76)   342     66        2 (1238/1458-1309/1493) -11.3 W11-W17 N31-N33
      2 1417 1402 (N28W25 L= 86)   478     94        3 (1383/1370-1451/1438)  53.1 W22-W27 N26-N31
      3  225 1277 (N11E63 L=359)  1320    332        2 ( 188/1180- 263/1395) -45.0 E69-E56 N05-N17
      4  430 1293 (N14E43 L= 19)   766    123        2 ( 407/1244- 464/1349) -54.5 E45-E39 N10-N17
      5 1222  782 (S12W14 L= 76)   563    224        1 (1195/ 693-1266/ 887) -61.4 W11-W18 S16-S05
      6  927  744 (S16E04 L= 57)   193    103        0 ( 890/ 710- 966/ 779)  45.0 E06-E01 S18-S13
      7 1674  641 (S18W50 L=112)   337     59        0 (1650/ 610-1700/ 663)  53.1 W48-W52 S20-S16
      8 1542  448 (S32W45 L=107)   565     99        2 (1500/ 423-1580/ 482)  18.4 W41-W48 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