Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1247 1803 (N45W53 L= 76)  1381    145        3 (1171/1761-1308/1835) -11.3 W45-W61 N40-N48
      2 1245 1440 (N24W26 L= 49)   373     67        1 (1213/1416-1267/1477)  21.8 W23-W28 N21-N25
      3 1059 1417 (N28W14 L= 37)   359     88        0 (1034/1378-1098/1451)  18.4 W11-W17 N25-N29
      4 1773 1299 (S00W58 L= 81)   371     40        0 (1750/1278-1800/1317) -11.3 W54-W60 S01-N00
      5  983 1197 (N17W02 L= 25)   831    229        2 ( 936/1104-1065/1290)  32.5 E02-W09 N13-N20
      6 1368 1220 (N08W25 L= 48)   220     53        0 (1352/1192-1386/1252)  45.0 W23-W26 N06-N09
      7  618 1127 (N23E22 L=  1)   536    187        0 ( 526/1086- 686/1155) -31.0 E28-E18 N18-N25
      8  372 1103 (N27E42 L=342)   379    114        1 ( 354/1045- 388/1149)  80.5 E43-E40 N23-N29
      9 1711  927 (S19W41 L= 64)   834    220        1 (1611/ 899-1816/ 953) -20.2 W32-W51 S22-S14
     10 1318  869 (S10W12 L= 35)   355    104        0 (1285/ 821-1342/ 917)  33.7 W09-W15 S11-S07
     11  458  833 (N09E40 L=343)   387    101        0 ( 433/ 790- 492/ 879) -80.5 E41-E38 N06-N12
     12  916  756 (S05E13 L= 10)   267    105        0 ( 889/ 710- 950/ 796)  38.7 E15-E10 S07-S03
     13 1197  629 (S20E01 L= 22)   304     99        2 (1178/ 591-1218/ 676)  78.7 E02-W01 S22-S17
     14  832  552 (S15E24 L=359)   615    164        2 ( 800/ 474- 876/ 617)  37.9 E28-E19 S19-S12
     15  649  394 (S21E44 L=339)   588    186        1 ( 605/ 354- 719/ 458) -83.7 E47-E39 S25-S16
     16  914  376 (S29E27 L=356)   641     94        1 ( 881/ 338- 953/ 409)  23.2 E30-E23 S31-S27
     17  574  349 (S23E53 L=330)   990     97        3 ( 531/ 323- 626/ 387) -45.0 E57-E49 S25-S19
  

    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