Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1750 (N58W00 L= 62)   366     37        1 ( 836/1726- 889/1767)   0.0 E02-W03 N56-N59
      2  708 1720 (N58E19 L= 43)   823    159        2 ( 637/1689- 778/1749)   6.3 E27-E09 N55-N60
      3 1025 1433 (N31W06 L= 68)   406    127        0 ( 967/1387-1059/1459) -45.0 W02-W08 N27-N33
      4  821 1358 (N29E09 L= 53)   345    126        0 ( 775/1326- 882/1399) -40.6 E11-E04 N26-N32
      5  234 1190 (N24E59 L=  3)   674    179        0 ( 200/1148- 280/1269) -33.7 E65-E52 N21-N29
      6  381 1175 (N22E43 L= 19)   644    181        1 ( 361/1096- 398/1257) -73.3 E45-E41 N17-N27
      7 1533 1232 (N11W36 L= 97)   332     69        1 (1499/1209-1567/1257) -26.6 W33-W37 N08-N12
      8 1591 1167 (N06W38 L=100)   298     90        1 (1575/1128-1613/1207)  78.7 W36-W40 N03-N08
      9 1901  962 (S13W67 L=128)  3030    222        4 (1865/ 873-1934/1059) -85.2 W61-W74 S19-S07
     10 1188  743 (S13W06 L= 68)   357    121        0 (1144/ 717-1240/ 782)  18.4 W03-W09 S14-S11
     11 1751  624 (S32W51 L=113)  1360    200        0 (1671/ 555-1833/ 676)  15.5 W42-W62 S35-S29
     12  382  634 (S11E50 L= 11)   299     78        1 ( 370/ 596- 398/ 664)  78.7 E51-E48 S13-S08
     13  234  607 (S12E70 L=352)   686     56        1 ( 205/ 578- 257/ 624) -31.0 E75-E65 S14-S11
     14 1261  359 (S42W07 L= 69)   434     75        0 (1229/ 341-1304/ 378)   0.0 W04-W10 S43-S39
     15 1068  266 (S46E11 L= 51)  1042    137        2 (1012/ 234-1141/ 294)  14.0 E16-E04 S48-S43
     16  899  202 (S50E30 L= 32)  1378    128        2 ( 856/ 157- 942/ 258)  32.0 E40-E23 S55-S44
  

    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