Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  690 1475 (N31E26 L= 39)   342     64        1 ( 653/1458- 729/1490)   0.0 E28-E22 N29-N31
      2 1743 1332 (N23W55 L=120)   361    118        1 (1692/1299-1790/1349) -18.4 W50-W60 N20-N23
      3 1802 1276 (N19W61 L=125)   380     72        2 (1781/1236-1818/1306)  33.7 W56-W63 N16-N20
      4  927 1252 (N17E07 L= 58)   402    125        1 ( 870/1226- 980/1272)   8.1 E10-E03 N15-N18
      5 1672 1210 (N15W45 L=110)   291     87        1 (1628/1189-1704/1228)   0.0 W41-W48 N13-N16
      6  155  820 (S14E72 L=352)  1554    119        4 ( 123/ 763- 185/ 869) -45.0 E79-E68 S17-S10
      7  274  762 (S16E56 L=  9)   557    116        2 ( 248/ 711- 297/ 806)  50.2 E60-E52 S19-S13
      8 1791  734 (S15W59 L=124)   913    156        1 (1759/ 670-1829/ 796)  76.0 W56-W62 S19-S11
      9  584  619 (S24E30 L= 35)   469     93        0 ( 545/ 582- 617/ 655) -38.7 E32-E27 S26-S21
     10 1273  448 (S34W20 L= 85)   943    233        4 (1236/ 342-1313/ 550) -67.6 W17-W25 S43-S26
  

    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