Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1690 (N37E41 L=355)  2067    280        3 ( 477/1598- 704/1741) -17.8 E54-E26 N30-N41
      2 1130 1667 (N38W03 L= 39)   322     68        0 (1091/1653-1170/1683) -18.4 W00-W06 N37-N39
      3 1079 1561 (N30W00 L= 36)   271     99        0 (1031/1539-1126/1587) -18.4 E03-W03 N28-N31
      4  883 1396 (N17E11 L= 25)   711    200        0 ( 829/1312- 927/1481)  76.0 E13-E08 N11-N23
      5 1804 1412 (N26W60 L= 95)  1178     83        1 (1764/1373-1834/1447) -23.2 W54-W64 N24-N29
      6  318 1375 (N14E53 L=343)   994    206        2 ( 250/1321- 390/1437) -21.3 E62-E44 N10-N18
      7  258 1310 (N10E57 L=339)  1209    202        2 ( 200/1235- 289/1378) -34.7 E65-E52 N05-N14
      8 1077 1306 (N13W02 L= 37)   554    234        0 ( 976/1235-1161/1359) -24.4 E04-W07 N08-N15
      9  171 1086 (S03E63 L=332)  1334    293        1 ( 141/ 972- 219/1219)  65.2 E67-E58 S09-N04
     10  366 1097 (S03E44 L=352)   414    104        0 ( 332/1055- 398/1140)  45.0 E46-E41 S05-N00
     11 1838  961 (S02W58 L= 94)  1392    160        2 (1809/ 881-1870/1035) -49.4 W54-W63 S05-N02
     12 1110  990 (S06W05 L= 41)   208     38        0 (1084/ 972-1133/1009)  33.7 W03-W06 S07-S04
     13  820  859 (S16E12 L= 24)   343    146        1 ( 745/ 837- 882/ 880) -12.5 E16-E07 S16-S14
     14  574  820 (S19E28 L=  7)   226     81        2 ( 542/ 799- 617/ 847) -21.8 E30-E25 S20-S17
     15 1659  695 (S19W47 L= 83)   950     49        0 (1608/ 653-1704/ 733)  35.5 W42-W50 S21-S16
     16 1700  658 (S21W52 L= 88)  1051    228        4 (1618/ 577-1781/ 723)  45.0 W46-W58 S26-S16
     17 1479  685 (S22W33 L= 69)   258     93        0 (1438/ 658-1513/ 717) -15.9 W29-W36 S23-S20
     18 1467  390 (S41W44 L= 80)   573     59        1 (1439/ 353-1495/ 413) -29.7 W40-W47 S43-S39
  

    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