Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1772 (N37E34 L= 16)  2817    583        4 ( 863/1649- 971/1925) -32.3 E53-E23 N28-N49
      2  402 1515 (N07E56 L=354)   526    126        0 ( 388/1462- 419/1565) -50.2 E59-E53 N04-N10
      3 1610 1509 (N43W26 L= 77)   331     60        1 (1591/1480-1632/1527)  53.1 W24-W28 N40-N44
      4 1358 1407 (N26W09 L= 59)   353    118        0 (1327/1343-1376/1448) -51.3 W06-W10 N21-N28
      5 1806 1370 (N41W47 L= 97)  1633    219        2 (1738/1308-1863/1423) -14.0 W39-W54 N38-N43
      6  785 1340 (N06E22 L= 28)   412     44        0 ( 761/1306- 810/1374) -63.4 E23-E20 N04-N08
      7 1852 1257 (N34W52 L=102)   756    148        1 (1808/1231-1877/1286)  36.9 W47-W55 N30-N37
      8 1137 1138 (N04W03 L= 53)   203     26        0 (1119/1118-1154/1162) -71.6 W02-W04 N02-N05
      9 1049 1150 (N02E02 L= 48)   205    107        0 ( 995/1140-1106/1159)  31.0 E04-W01 N00-N03
     10 1673 1033 (N14W39 L= 89)   864    217        0 (1635/ 959-1728/1081)  76.0 W36-W43 N09-N17
     11 1586 1024 (N11W33 L= 83)   470     84        1 (1560/ 979-1615/1070) -80.5 W30-W35 N07-N13
     12  842  941 (S15E08 L= 42)   305    108        1 ( 819/ 886- 867/ 996)  90.0 E09-E07 S18-S10
     13 1697  890 (N07W45 L= 95)   807    140        1 (1644/ 837-1736/ 946)  66.0 W41-W48 N02-N11
     14 1835  892 (N12W58 L=108)   814    117        2 (1809/ 832-1862/ 939)  82.9 W55-W61 N07-N15
     15  386  894 (S29E38 L= 12)   727    162        1 ( 321/ 845- 445/ 935)  38.7 E43-E33 S33-S25
     16  690  850 (S24E15 L= 35)   583    181        0 ( 611/ 819- 785/ 896)  45.0 E19-E09 S27-S18
     17  868  699 (S28W00 L= 50)   456     74        2 ( 834/ 666- 901/ 732)  63.4 E01-W02 S29-S25
     18  649  625 (S39E13 L= 37)   433     70        0 ( 615/ 607- 693/ 650)  18.4 E15-E09 S40-S36
     19 1186  296 (S42W40 L= 90)   367     43        0 (1157/ 276-1214/ 316)  11.3 W37-W42 S43-S41
  

    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