Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1007 1754 (N51W30 L= 10)   300     37        1 ( 983/1731-1033/1769) -11.3 W26-W32 N49-N52
      2 1533 1644 (N25W58 L= 38)   765     73        1 (1496/1613-1558/1673) -78.7 W55-W61 N22-N27
      3 1538 1533 (N19W50 L= 30)   347     70        1 (1519/1498-1555/1568)  26.6 W47-W53 N17-N20
      4 1301 1510 (N25W32 L= 13)   401    111        0 (1245/1489-1354/1533) -39.8 W29-W35 N22-N28
      5 1202 1401 (N22W22 L=  2)   650    225        0 (1103/1351-1283/1437) -48.4 W16-W24 N17-N26
      6  244 1405 (N49E62 L=278)   512     58        0 ( 218/1379- 274/1422)   6.3 E67-E56 N47-N50
      7  972 1216 (N19W02 L=342)   936    361        1 ( 884/1084-1096/1335) -75.3 W00-W05 N08-N27
      8  410 1259 (N36E36 L=304)  1174    198        2 ( 338/1208- 501/1312) -29.4 E44-E28 N31-N40
      9  579 1135 (N24E24 L=316)   336     77        0 ( 548/1106- 619/1162) -45.0 E26-E21 N21-N26
     10  174 1021 (N25E65 L=275)   667     66        1 ( 156/ 989- 191/1060) -63.4 E67-E62 N23-N27
     11  204  930 (N20E62 L=279)  2790    212        8 ( 137/ 853- 251/1039)  51.3 E71-E55 N14-N25
     12 1466  886 (S13W22 L=  2)   952    302        0 (1405/ 770-1523/1028)  42.7 W15-W28 S18-S06
     13 1746  936 (S20W44 L= 24)   906    172        2 (1725/ 850-1762/1025)  63.4 W41-W46 S24-S14
     14  707  914 (N08E21 L=319)   613    169        0 ( 643/ 858- 771/ 979) -66.0 E23-E18 N04-N13
     15  412  888 (N13E42 L=298)   293     66        1 ( 376/ 865- 443/ 908) -31.0 E44-E39 N11-N14
     16 1206  635 (S20E00 L=340)   295     99        1 (1182/ 611-1252/ 671) -63.4 E01-W02 S22-S17
     17 1466  615 (S30W16 L=356)   551    104        1 (1440/ 569-1506/ 648)  33.7 W13-W19 S32-S27
     18  981  454 (S26E19 L=321)   449    171        0 ( 933/ 395-1043/ 516)  20.0 E24-E13 S28-S23
     19  958  267 (S38E31 L=309)   726    124        0 ( 930/ 208- 982/ 326)  51.3 E36-E27 S43-S33
  

    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