Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  985 1864 (N48E38 L=360)   516     88        2 ( 960/1844-1008/1888) -56.3 E40-E34 N44-N50
      2 1020 1773 (N41E27 L= 10)  1229    290        4 ( 984/1696-1049/1850) -49.8 E33-E21 N34-N47
      3 1077 1638 (N33E16 L= 21)   810    145        2 (1044/1589-1132/1696) -18.4 E20-E11 N29-N36
      4  377 1573 (N11E62 L=336)  3533    521       12 ( 269/1464- 486/1685)  64.8 E70-E51 N02-N19
      5  521 1504 (N09E45 L=352)   591    190        1 ( 488/1430- 543/1567) -37.9 E50-E40 N05-N12
      6 1745 1428 (N42W41 L= 78)   902    216        1 (1674/1373-1811/1480)  -7.6 W33-W48 N39-N44
      7 1516 1324 (N26W21 L= 59)   255     88        0 (1500/1272-1530/1359) -80.5 W19-W22 N22-N28
      8  979 1271 (N07E09 L= 28)   333     89        0 ( 956/1226-1003/1325) -80.5 E10-E08 N04-N10
      9  334 1223 (S11E50 L=348)   311     48        0 ( 313/1199- 356/1249)  63.4 E51-E48 S12-S08
     10 1238 1056 (N02W11 L= 48)   276    123        0 (1182/1041-1294/1084)  38.7 W07-W14 N00-N05
     11 1326 1043 (N04W16 L= 54)   203     22        0 (1308/1022-1342/1066)  71.6 W15-W17 N02-N05
     12 1790  904 (N11W53 L= 90)  3250    669        6 (1689/ 737-1859/1012) -11.7 W41-W66 N01-N18
     13 1889  808 (N10W69 L=106)   956    109        1 (1874/ 761-1909/ 854) -80.5 W66-W72 N06-N12
     14  855  761 (S24E03 L= 35)   608    186        0 ( 768/ 735- 952/ 803)  41.6 E07-W02 S27-S19
     15 1015  606 (S29W12 L= 50)   289     56        2 ( 982/ 585-1047/ 633)  14.0 W10-W14 S29-S27
     16  789  544 (S40E00 L= 37)   491    134        0 ( 720/ 529- 849/ 558)  26.6 E05-W03 S41-S37
     17 1292  234 (S42W52 L= 89)   560     63        0 (1256/ 218-1325/ 247)  45.0 W48-W54 S43-S40
  

    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