Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  331 1495 (N07E59 L=137)   623     92        3 ( 303/1452- 358/1532) -68.2 E62-E56 N04-N09
      2 1537 1364 (N25W23 L=220)   302     70        1 (1503/1339-1569/1387)  45.0 W21-W25 N22-N26
      3  404 1280 (S05E45 L=152)   252     48        1 ( 386/1254- 424/1311)  90.0 E45-E42 S06-S03
      4 1152 1184 (N05W03 L=200)   238     73        0 (1115/1158-1195/1216)  51.3 W01-W05 N02-N07
      5  669 1110 (S10E23 L=174)   276    106        0 ( 647/1049- 682/1155) -63.4 E24-E21 S13-S07
      6  528 1008 (S18E30 L=167)   223     96        0 ( 488/ 977- 558/1045)  51.3 E32-E28 S20-S15
      7 1882  870 (N08W64 L=260)   444     95        2 (1871/ 837-1898/ 919) -68.2 W61-W66 N05-N10
      8 1480  566 (S21W41 L=238)   429    103        2 (1445/ 527-1506/ 605) -78.7 W38-W42 S23-S18
      9 1281  567 (S26W28 L=225)   283     52        2 (1250/ 547-1318/ 585)  26.6 W26-W30 S27-S24
     10  862  215 (S65W22 L=218)   667    127        3 ( 831/ 189- 889/ 248)  90.0 W17-W26 S67-S60
  

    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