Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1100 1716 (N38E19 L= 67)   290     26        0 (1080/1695-1123/1733) -18.4 E20-E17 N36-N39
      2 1285 1493 (N29W02 L= 87)  1170    251        1 (1197/1395-1331/1587) -36.9 E06-W06 N23-N33
      3 1071 1329 (N12E06 L= 80)   321    106        1 (1040/1277-1090/1375) -45.0 E08-E03 N09-N14
      4 1128 1248 (N09E00 L= 85)   225     81        2 (1106/1211-1165/1273)  90.0 E02-W01 N06-N11
      5 1613 1240 (N24W30 L=115)   395     75        0 (1573/1213-1641/1263)   9.5 W26-W32 N22-N25
      6 1308 1185 (N11W11 L= 97)   721    266        0 (1203/1094-1352/1255) -17.1 W03-W16 N06-N13
      7 1590 1093 (N14W31 L=117)   405     83        0 (1572/1052-1607/1139) -80.5 W30-W32 N11-N17
      8 1689  980 (N11W42 L=127)   661    194        0 (1670/ 900-1720/1062) -60.9 W39-W45 N07-N16
      9 1791  997 (N16W50 L=136)  1070    166        1 (1737/ 932-1846/1051) -12.1 W43-W58 N14-N17
     10 1211  518 (S28W28 L=114)   392     79        0 (1172/ 502-1256/ 532)  36.9 W25-W30 S29-S26
     11 1666  447 (S16W67 L=152)  1830    180        7 (1625/ 370-1708/ 520) -51.3 W61-W74 S21-S11
     12  754  497 (S45E00 L= 85)   295     38        0 ( 735/ 474- 781/ 519)  76.0 E02-W01 S46-S42
     13  625  480 (S51E11 L= 75)   519     76        0 ( 596/ 454- 659/ 513) -14.0 E14-E06 S52-S48
     14 1559  401 (S23W60 L=146)   436     91        1 (1538/ 366-1593/ 430) -78.7 W57-W62 S25-S20
     15 1197  373 (S37W36 L=122)   446    138        1 (1145/ 322-1250/ 409)  90.0 W34-W37 S41-S33
     16 1321  342 (S34W47 L=133)   993    292        0 (1194/ 284-1446/ 389)  59.0 W43-W52 S42-S27
     17 1449  300 (S32W61 L=147)  1244    169        2 (1411/ 242-1486/ 365) -14.7 W54-W74 S34-S27
  

    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