Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1251 1806 (N45W53 L= 75)  1169    158        3 (1165/1765-1316/1827) -19.7 W46-W61 N40-N48
      2 1260 1448 (N24W27 L= 49)   378     73        2 (1227/1422-1287/1480)   9.5 W24-W30 N21-N25
      3 1066 1416 (N27W14 L= 36)   275     42        0 (1048/1385-1087/1443)  36.9 W12-W16 N25-N28
      4 1784 1306 (S01W59 L= 81)   444     53        0 (1761/1281-1810/1325)  -9.5 W56-W62 S01-N00
      5 1004 1213 (N18W04 L= 26)   763    259        1 ( 958/1110-1084/1313)  33.7 E01-W11 N13-N21
      6 1379 1232 (N08W26 L= 48)   285     78        1 (1354/1201-1402/1266) -90.0 W24-W27 N06-N10
      7  632 1131 (N23E21 L=  1)   598    209        0 ( 533/1086- 702/1159) -35.0 E27-E17 N18-N25
      8 1719  934 (S18W41 L= 63)   820    213        1 (1625/ 908-1824/ 958) -24.0 W34-W52 S22-S14
      9 1334  879 (S10W14 L= 36)   355     98        0 (1303/ 833-1361/ 923)  33.7 W10-W16 S11-S07
     10  465  842 (N10E39 L=343)   315     87        0 ( 441/ 809- 493/ 890) -68.2 E40-E37 N07-N12
     11  204  785 (N11E66 L=316)   468     99        1 ( 181/ 739- 233/ 824) -71.6 E68-E63 N07-N13
     12  931  765 (S05E12 L= 10)   253    100        0 ( 909/ 722- 968/ 799)  31.0 E14-E09 S07-S03
     13 1210  637 (S20W00 L= 22)   352     91        3 (1192/ 595-1237/ 690)  80.5 E01-W01 S22-S16
     14  847  552 (S16E23 L=359)   697    223        0 ( 808/ 472- 905/ 634)  32.5 E28-E17 S20-S11
     15  629  415 (S19E44 L=338)   387     98        0 ( 614/ 371- 648/ 462)  50.2 E46-E40 S22-S16
     16  923  375 (S29E26 L=356)   645    103        1 ( 892/ 329- 970/ 403)  20.6 E30-E22 S31-S27
     17  581  356 (S22E52 L=330)   572     50        1 ( 562/ 329- 601/ 388)  59.0 E54-E48 S24-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