Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  913 1618 (N45W13 L= 74)   782    134        0 ( 855/1585- 979/1641) -32.0 W08-W17 N41-N48
      2 1588 1488 (N16W51 L=113)   889    240        3 (1559/1380-1626/1578)  32.7 W45-W59 N09-N19
      3 1460 1371 (N14W36 L= 98)   491    119        1 (1398/1347-1515/1394) -18.4 W31-W40 N12-N15
      4  849 1277 (N26E04 L= 58)   226    114        0 ( 790/1266- 899/1291) -15.9 E07-W00 N24-N26
      5  782 1196 (N23E10 L= 52)   322     63        1 ( 750/1164- 804/1226)  63.4 E11-E07 N20-N24
      6 1755 1107 (S10W49 L=110)   479    107        1 (1712/1066-1792/1144)   5.7 W43-W53 S10-S08
      7  471  983 (N18E35 L= 26)   607    173        0 ( 415/ 935- 535/1041)  20.0 E40-E29 N15-N19
      8  187 1009 (N24E63 L=358)   500     84        0 ( 150/ 986- 221/1028)   0.0 E69-E58 N23-N25
      9  254  940 (N19E56 L=  6)   590     91        2 ( 226/ 908- 295/ 984) -50.2 E58-E51 N16-N22
     10  786  931 (N08E16 L= 46)   288     78        0 ( 769/ 884- 804/ 968)  53.1 E17-E14 N05-N09
     11  158  873 (N16E70 L=351)   903     84        2 ( 140/ 831- 175/ 920) -80.5 E73-E67 N13-N19
     12 1323  834 (S12W12 L= 73)   351    119        0 (1300/ 795-1348/ 885)  68.2 W10-W13 S13-S08
     13 1605  782 (S23W30 L= 91)   607    157        0 (1560/ 716-1647/ 836)  32.0 W25-W33 S26-S21
     14 1690  740 (S29W36 L= 98)   413     86        0 (1662/ 709-1713/ 790)  38.7 W33-W38 S30-S26
     15 1175  659 (S18E01 L= 61)  1059    416        4 ( 995/ 633-1371/ 693) -20.0 E11-W11 S22-S13
     16 1711  591 (S41W39 L=100)   539     54        1 (1685/ 566-1738/ 624)  45.0 W35-W42 S42-S38
     17  406  426 (S15E63 L=358)  1284    159        4 ( 383/ 361- 425/ 505)  29.1 E74-E56 S19-S09
     18  690  423 (S20E39 L= 23)   365     61        1 ( 673/ 393- 713/ 457)  78.7 E40-E36 S22-S17
  

    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