Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  434 1442 (N30E38 L=193)   944    247        1 ( 350/1394- 524/1489) -17.4 E46-E30 N26-N34
      2 1317 1454 (N22W23 L=255)   714    202        0 (1226/1431-1419/1478)   4.1 W16-W30 N20-N23
      3 1020 1342 (N18W03 L=234)   259     74        0 ( 975/1326-1064/1358) -21.8 W00-W05 N16-N18
      4  298 1290 (N22E48 L=183)   347     74        0 ( 273/1264- 323/1313) -38.7 E51-E46 N19-N23
      5 1909 1278 (N07W73 L=304)   601    115        2 (1881/1261-1927/1296)   0.0 W67-W77 N05-N08
      6  266 1182 (N15E51 L=181)   346     62        1 ( 247/1152- 284/1215) -76.0 E52-E49 N13-N17
      7 1237 1139 (N04W14 L=245)   334     68        0 (1220/1095-1253/1182) -90.0 W12-W14 N01-N06
      8  300 1084 (N09E48 L=184)   671     95        2 ( 265/1042- 332/1112) -38.7 E51-E45 N06-N10
      9 1558 1078 (S02W34 L=266)   626    102        1 (1505/1039-1619/1110) -20.6 W30-W38 S05-N00
     10  157  885 (S01E66 L=165)  1320    333        2 ( 101/ 789- 209/ 948)  21.3 E79-E58 S06-N02
     11 1385  884 (S13W21 L=253)   352     66        1 (1359/ 844-1407/ 915) -90.0 W19-W22 S15-S10
     12  867  863 (S09E11 L=221)   239     31        0 ( 838/ 844- 890/ 881)   0.0 E12-E09 S10-S08
     13  217  709 (S12E63 L=169)  3320    469       10 ( 147/ 609- 277/ 837) -75.1 E69-E56 S18-S03
     14 1170  796 (S16W07 L=239)   192     18        0 (1150/ 780-1191/ 815) -33.7 W05-W08 S17-S14
     15 1262  739 (S20W12 L=244)   488    152        0 (1211/ 683-1301/ 811)  66.8 W08-W15 S23-S15
     16 1364  722 (S22W19 L=251)   632    183        0 (1293/ 637-1438/ 799)  42.0 W14-W24 S27-S18
     17  482  535 (S26E44 L=188)   852    178        1 ( 425/ 498- 541/ 583)  -5.7 E48-E38 S28-S23
     18 1391  388 (S47W24 L=255)  1387    236        3 (1293/ 349-1482/ 446)  14.0 W15-W31 S49-S42
  

    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