Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1160 1550 (N37W05 L=304)   671    217        1 (1057/1495-1271/1593)  -3.8 E02-W13 N33-N39
      2 1699 1485 (N36W51 L=349)   872    189        2 (1633/1439-1771/1516) -22.6 W44-W57 N33-N38
      3  529 1449 (N24E39 L=259)   300     31        0 ( 509/1426- 547/1474)  71.6 E40-E37 N22-N25
      4 1410 1430 (N31W24 L=322)   422     95        1 (1375/1390-1446/1465) -78.7 W21-W27 N27-N32
      5  634 1440 (N25E31 L=267)   528    130        0 ( 579/1411- 687/1464) -12.5 E35-E26 N23-N25
      6 1594 1319 (N24W38 L=336)  1668    295        5 (1482/1238-1698/1385) -25.2 W29-W46 N19-N28
      7  231  983 (S08E57 L=241)  1439    224        1 ( 201/ 876- 255/1096) -74.1 E62-E54 S14-N00
      8  736  735 (S19E16 L=282)   284     51        1 ( 704/ 710- 766/ 757)  36.9 E18-E14 S20-S16
      9 1787  660 (S16W62 L=360)   792    135        2 (1751/ 631-1812/ 702)  31.0 W57-W64 S17-S12
  

    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