Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1432 1504 (N21W40 L= 15)   619    119        2 (1400/1451-1469/1555)  24.0 W36-W45 N18-N23
      2  905 1493 (N37W07 L=341)   293     91        0 ( 864/1465- 942/1521)   0.0 W03-W10 N35-N37
      3 1196 1485 (N28W24 L=359)   294    107        0 (1149/1465-1253/1502) -15.9 W20-W28 N26-N28
      4  718 1450 (N40E08 L=327)   538    126        1 ( 672/1419- 772/1472) -35.5 E11-E04 N36-N41
      5 1650 1294 (N04W46 L= 20)   412     76        1 (1614/1273-1685/1325)   0.0 W42-W48 N02-N05
      6  169 1264 (N40E68 L=266)   823     77        0 ( 126/1242- 206/1286)  -8.1 E81-E60 N38-N41
      7  634 1188 (N26E19 L=315)  1482    304        3 ( 485/1146- 777/1211) -24.2 E30-E10 N20-N30
      8  221  869 (N16E60 L=275)  4550    668        7 ( 120/ 664- 470/1008)  19.9 E75-E35 N03-N23
      9 1659  922 (S16W36 L= 10)   867    215        2 (1615/ 859-1695/1002)  36.9 W31-W40 S19-S12
     10 1427  907 (S10W20 L=354)   853    187        0 (1369/ 833-1469/ 988)  33.7 W14-W23 S13-S06
     11 1173  669 (S17E01 L=334)   350    107        0 (1144/ 622-1205/ 721)  38.7 E03-W02 S18-S14
     12  885  643 (S11E18 L=317)   274    107        2 ( 850/ 589- 906/ 684) -90.0 E19-E15 S14-S08
     13  950  657 (S12E14 L=321)   335     66        2 ( 917/ 629- 983/ 683)  26.6 E15-E11 S12-S10
     14  986  620 (S15E13 L=322)   247     12        0 ( 964/ 600-1011/ 640) -18.4 E14-E11 S15-S13
     15 1137  477 (S28E09 L=326)   353    111        1 (1124/ 414-1156/ 528)  54.5 E11-E05 S31-S24
     16  407  386 (S17E65 L=269)  3110    203        7 ( 344/ 289- 455/ 454)  41.2 E84-E58 S25-S11
     17  615  368 (S22E47 L=288)   866    112        0 ( 589/ 303- 650/ 415)  28.6 E53-E41 S25-S18
     18 1192  314 (S42E12 L=323)   585    102        0 (1174/ 270-1208/ 378)  52.1 E15-E08 S45-S36
     19 1409  266 (S56W03 L=338)  1434    160        2 (1335/ 248-1487/ 291) -21.8 E04-W12 S61-S50
  

    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