Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1092 1547 (N30W08 L=170)   687    171        0 (1020/1489-1159/1596) -45.0 W03-W11 N25-N33
      2  472 1484 (N29E37 L=125)   868    156        1 ( 416/1444- 544/1536) -28.3 E43-E30 N26-N33
      3  270 1305 (N19E53 L=109)  3825    386        9 ( 178/1164- 345/1465)  68.7 E61-E44 N11-N29
      4 1590 1327 (N13W38 L=200)   283    100        0 (1542/1312-1640/1353) -15.9 W35-W42 N12-N15
      5  198 1275 (N18E60 L=102)   688    120        5 ( 175/1228- 234/1330)  90.0 E63-E55 N15-N21
      6 1488 1209 (N07W30 L=192)   233     68        0 (1451/1197-1527/1222) -11.3 W27-W32 N05-N07
      7 1267 1100 (N01W15 L=177)   711    320        1 (1141/1049-1403/1129)  -3.6 W07-W23 S02-N02
      8  455  851 (S09E37 L=125)  1064    292        2 ( 403/ 718- 521/ 932)  47.5 E43-E31 S16-S04
      9 1224  873 (S12W11 L=173)   484    124        1 (1179/ 840-1293/ 904) -15.9 W08-W15 S13-S10
     10 1030  823 (S14E01 L=161)  1244    383        2 ( 905/ 758-1157/ 879) -19.7 E08-W06 S18-S10
     11  224  792 (S11E59 L=103)   967    102        0 ( 186/ 750- 261/ 839)  20.0 E65-E54 S13-S08
     12  511  771 (S14E34 L=128)   285    102        1 ( 484/ 721- 534/ 807)  45.0 E37-E32 S17-S12
     13 1090  640 (S26W02 L=164)   217     61        1 (1063/ 615-1118/ 673)  36.9 W00-W04 S27-S23
     14 1469  526 (S35W30 L=192)   319     84        0 (1429/ 501-1501/ 560)  36.9 W27-W32 S36-S33
     15  396  503 (S30E53 L=109)   549     53        0 ( 369/ 484- 429/ 528) -45.0 E56-E50 S31-S28
     16  345  438 (S34E65 L= 97)  1057    132        2 ( 325/ 371- 372/ 502)  26.6 E76-E57 S38-S29
     17 1665  469 (S40W53 L=215)  1015    115        1 (1600/ 445-1722/ 490)   0.0 W45-W62 S41-S38
  

    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