Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1653 (N40W30 L= 33)   266     85        1 (1101/1625-1170/1674) -59.0 W28-W31 N37-N42
      2  416 1317 (N39E36 L=327)   387     97        1 ( 373/1285- 456/1348) -29.1 E40-E31 N36-N41
      3  902 1307 (N26W00 L=  3)   210     65        1 ( 867/1293- 937/1322) -18.4 E01-W02 N25-N27
      4 1853 1269 (S04W64 L= 67)   420     39        0 (1826/1247-1873/1286) -31.0 W61-W66 S05-S02
      5 1387 1220 (N08W26 L= 28)   359     78        0 (1357/1189-1417/1252)  21.8 W23-W28 N06-N09
      6 1330 1130 (N05W20 L= 22)   188     39        0 (1314/1108-1346/1156)  71.6 W18-W20 N03-N06
      7  304 1068 (N26E48 L=315)  1726    254        4 ( 198/1003- 414/1125) -25.5 E60-E38 N19-N30
      8  148  903 (N18E70 L=293)  1381    183        1 (  87/ 852- 210/ 951) -15.3 E83-E61 N14-N21
      9 1407  862 (S12W18 L= 20)   478    274        1 (1312/ 832-1521/ 891) -12.1 W11-W25 S15-S08
     10  975  685 (S11E11 L=351)   641    211        1 ( 921/ 614-1021/ 785)  41.6 E16-E06 S13-S05
     11  534  505 (S11E45 L=318)   330    140        0 ( 492/ 474- 576/ 578) -90.0 E46-E42 S13-S06
     12  633  545 (S11E36 L=326)   285     74        1 ( 594/ 531- 672/ 561)  -9.5 E39-E33 S12-S09
     13  754  490 (S17E31 L=332)   473    129        2 ( 724/ 434- 789/ 552)  41.2 E34-E26 S20-S13
     14 1583  509 (S41W24 L= 27)   370     65        1 (1568/ 475-1604/ 537) -90.0 W22-W26 S43-S38
     15 1462  417 (S43W12 L= 15)   532    125        0 (1414/ 378-1510/ 472)  16.7 W07-W17 S45-S40
     16  790  316 (S29E38 L=325)   332     98        2 ( 769/ 267- 803/ 355)  49.4 E41-E35 S33-S26
     17  874  237 (S38E38 L=324)  1208    187        0 ( 834/ 152- 907/ 317)  53.1 E46-E33 S46-S30
     18 1210  173 (S57E21 L=342)  1768    196        3 (1130/ 134-1304/ 219)  -5.7 E31-E09 S59-S52
  

    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