Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1848 (N63W32 L= 72)   748     95        1 ( 861/1809- 922/1888)  13.2 W25-W42 N59-N65
      2  855 1785 (N59W20 L= 59)   339     53        1 ( 823/1767- 883/1805) -31.0 W16-W22 N56-N60
      3 1242 1729 (N40W42 L= 82)   297     48        0 (1218/1706-1269/1753) -90.0 W40-W44 N37-N41
      4  698 1722 (N60W00 L= 39)   456     84        1 ( 664/1692- 738/1746)   6.3 E04-W05 N58-N61
      5 1533 1713 (N28W63 L=103)   597     49        0 (1507/1689-1562/1737) -56.3 W59-W66 N26-N29
      6 1436 1548 (N23W43 L= 82)  1348    314        3 (1276/1503-1573/1595) -39.1 W34-W51 N17-N30
      7 1761 1391 (N05W60 L= 99)   794    161        1 (1742/1320-1781/1477)  45.0 W56-W65 N00-N08
      8 1237 1260 (N14W18 L= 58)   225    121        0 (1208/1213-1251/1309)  71.6 W16-W19 N11-N17
      9 1068 1237 (N17W08 L= 48)   855    262        4 ( 975/1185-1152/1289) -20.0 W02-W13 N14-N20
     10  515 1223 (N30E27 L= 12)   307     91        0 ( 468/1209- 561/1241)  -7.1 E31-E23 N29-N31
     11 1278 1166 (N08W18 L= 57)   279     82        0 (1237/1149-1302/1196)   0.0 W15-W19 N06-N08
     12  921 1157 (N16E02 L= 37)   408    161        0 ( 844/1142-1002/1169) -35.5 E06-W02 N13-N18
     13  419 1122 (N27E36 L=  3)   562    127        0 ( 361/1105- 497/1144) -20.0 E41-E30 N23-N28
     14 1181  897 (S04W05 L= 45)   437    114        0 (1134/ 877-1231/ 933)   9.5 W02-W08 S05-S03
     15  926  717 (S08E13 L= 26)   887    288        0 ( 881/ 604- 971/ 808)  77.0 E15-E09 S15-S02
     16  761  472 (S18E31 L=  8)  1150    341        1 ( 656/ 427- 888/ 532) -31.6 E36-E23 S23-S13
     17 1542  496 (S41W20 L= 59)   723     85        3 (1499/ 465-1583/ 521)   7.1 W15-W23 S42-S38
     18  475  359 (S20E60 L=339)  1805    256        1 ( 452/ 272- 502/ 478)  30.3 E73-E49 S26-S12
     19 1146  256 (S45E19 L= 21)  1791    332        1 (1025/ 173-1291/ 345)   8.1 E33-E04 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