Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1465 (N35W04 L=298)  1183    313        2 ( 791/1435-1077/1496) -17.5 E05-W14 N31-N38
      2 1524 1463 (N20W44 L=338)   319     67        0 (1506/1430-1541/1495)  68.2 W43-W45 N17-N22
      3 1275 1424 (N24W26 L=319)   814    291        1 (1164/1331-1368/1464) -50.7 W19-W29 N16-N28
      4 1260 1260 (N15W19 L=313)   584    197        0 (1219/1194-1317/1302)  56.3 W16-W22 N11-N18
      5 1785 1099 (S07W52 L=346)   410     96        1 (1759/1064-1802/1155)  51.3 W49-W54 S09-S04
      6 1837 1092 (S09W58 L=352)   394     53        1 (1811/1072-1868/1110)  -7.1 W54-W62 S10-S08
      7  618 1071 (N18E24 L=270)   319     56        1 ( 594/1042- 651/1101) -45.0 E25-E21 N15-N19
      8  211 1014 (N20E60 L=233)   876    135        1 ( 180/ 962- 263/1060) -37.9 E64-E54 N16-N23
      9 1931 1020 (S17W71 L=  5)   810     89        2 (1893/ 994-1953/1043)  -3.8 W63-W78 S19-S15
     10 1712  994 (S11W43 L=337)   355    100        1 (1681/ 968-1757/1019)   0.0 W40-W47 S12-S09
     11  142  969 (N18E71 L=222)   767     66        1 ( 121/ 933- 160/1001)  53.1 E75-E67 N15-N19
     12  639  887 (N07E26 L=268)   689    192        0 ( 601/ 798- 677/ 972) -85.2 E27-E24 N00-N12
     13 1674  755 (S25W37 L=331)   479    122        1 (1640/ 722-1709/ 810)  26.6 W34-W40 S26-S22
     14  271  613 (S05E65 L=229)  5624    464        8 ( 192/ 443- 365/ 775)  81.9 E82-E51 S17-N05
     15 1169  635 (S19W00 L=294)   228     80        0 (1129/ 620-1208/ 651) -14.0 E02-W02 S19-S17
     16 1192  544 (S26E00 L=293)   300     76        0 (1170/ 511-1213/ 581)  53.1 E02-W01 S27-S23
  

    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