Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  497 1719 (N22E62 L=  1)   861    155        3 ( 467/1660- 521/1799) -38.2 E72-E58 N17-N28
      2  436 1617 (N14E59 L=  4)  2105    186        7 ( 398/1541- 475/1706) -32.9 E71-E53 N09-N20
      3 1513 1478 (N36W18 L= 81)  2402    499        3 (1401/1255-1630/1665) -42.5 W06-W31 N25-N47
      4  516 1591 (N15E50 L= 13)   349     48        1 ( 489/1577- 551/1606)  26.6 E51-E47 N13-N15
      5  311 1506 (N05E64 L=359)  2260    256        4 ( 199/1460- 416/1567)  17.7 E81-E54 N01-N10
      6  526 1533 (N11E46 L= 17)   277     64        0 ( 491/1511- 552/1555)  53.1 E48-E43 N09-N13
      7  900 1439 (N15E18 L= 45)   391    126        0 ( 846/1399- 936/1472)   0.0 E21-E15 N12-N16
      8 1353 1391 (N25W09 L= 72)   376    140        0 (1295/1339-1408/1426)  -6.3 W04-W13 N23-N25
      9  635 1384 (N05E32 L= 31)   672    126        3 ( 576/1352- 683/1425)  50.2 E35-E29 N02-N08
     10 1782 1172 (N25W46 L=109)  1382    255        1 (1724/1093-1859/1272) -32.5 W39-W55 N21-N30
     11  247 1233 (S12E58 L=  5)   350     66        0 ( 232/1196- 258/1263) -36.9 E60-E56 S13-S10
     12 1466 1226 (N18W20 L= 83)   258     53        0 (1443/1196-1484/1254) -33.7 W18-W21 N16-N19
     13  551 1185 (S08E32 L= 31)   367    140        0 ( 519/1130- 573/1239) -56.3 E35-E30 S11-S05
     14  149 1022 (S26E67 L=356)  1705    133        5 ( 104/ 989- 212/1070)   9.0 E76-E57 S28-S22
     15 1152  868 (S10W11 L= 74)   284     56        1 (1117/ 842-1182/ 895)  53.1 W09-W12 S11-S07
     16 1873  763 (N06W68 L=131)  2319    176        6 (1828/ 699-1905/ 833) -63.4 W61-W77 N02-N10
     17 1723  702 (S02W54 L=117)  1938    199        3 (1689/ 610-1780/ 787) -28.3 W47-W60 S06-N01
  

    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