Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1086 1858 (N56W34 L=355)  1792    233        5 ( 985/1788-1173/1912)   4.4 W19-W47 N51-N60
      2 1487 1784 (N39W63 L= 24)  2512    317        4 (1354/1705-1612/1830) -31.0 W50-W77 N31-N45
      3 1687 1633 (N25W69 L= 30)   603     80        0 (1646/1607-1720/1658) -33.7 W64-W72 N22-N26
      4  945 1588 (N37W06 L=328)   954    224        2 ( 844/1541-1026/1612) -14.0 W00-W12 N32-N39
      5  524 1369 (N31E28 L=294)  1774    496        6 ( 312/1336- 777/1405) -17.1 E48-E09 N24-N37
      6 1215 1360 (N17W18 L=339)   615    189        2 (1161/1315-1284/1400)  -7.1 W14-W22 N14-N20
      7  291 1298 (N31E49 L=272)   391     77        0 ( 259/1259- 313/1324) -45.0 E53-E46 N27-N32
      8 1001 1277 (N16W03 L=325)   573    235        0 ( 891/1251-1105/1298) -22.6 E03-W09 N13-N19
      9  175 1249 (N29E64 L=257)   715    107        1 ( 144/1216- 211/1297) -23.2 E72-E57 N26-N32
     10  625 1232 (N20E21 L=300)   349    155        0 ( 553/1212- 701/1248) -24.0 E25-E16 N18-N22
     11  558 1128 (N15E27 L=295)   682    298        2 ( 451/1066- 674/1213) -42.9 E33-E19 N09-N22
     12  376  828 (N01E45 L=276)   379    101        0 ( 335/ 791- 419/ 874)  12.5 E49-E40 N00-N03
     13  546  644 (S13E37 L=284)   707    284        1 ( 451/ 581- 678/ 696) -39.8 E43-E30 S18-S08
     14  884  647 (S18E15 L=306)   228     43        1 ( 866/ 617- 903/ 676)  45.0 E17-E13 S19-S16
     15  726  440 (S28E33 L=289)   626    171        1 ( 705/ 361- 764/ 504)  48.4 E37-E27 S33-S24
     16 1234  457 (S38W03 L=324)   335    151        0 (1169/ 424-1298/ 498)  11.3 E02-W08 S39-S36
     17  889  458 (S30E21 L=301)   249    110        0 ( 850/ 428- 938/ 491)  14.0 E24-E16 S31-S29
     18  613  265 (S39E55 L=266)   869     86        2 ( 578/ 231- 644/ 306)  14.0 E64-E48 S40-S35
  

    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