Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  910 1546 (N28E10 L=224)   273    101        1 ( 889/1491- 936/1587) -66.8 E12-E08 N23-N30
      2 1157 1510 (N26W06 L=241)   377    136        1 (1081/1486-1219/1526)  -6.3 W01-W10 N25-N27
      3  788 1485 (N23E18 L=217)   670    168        1 ( 701/1446- 854/1508) -10.3 E24-E13 N20-N24
      4  694 1410 (N18E23 L=212)   227     83        0 ( 655/1384- 729/1430) -21.8 E25-E20 N16-N18
      5  343 1370 (N15E49 L=186)   308     80        0 ( 327/1330- 354/1406) -63.4 E50-E47 N12-N16
      6  581 1358 (N14E30 L=205)   290     57        0 ( 568/1320- 599/1390) -78.7 E31-E28 N11-N16
      7 1847 1267 (N17W59 L=294)  2655    541        6 (1713/1130-1966/1375) -16.7 W46-W76 N09-N22
      8 1192 1262 (N11W09 L=244)   399    108        0 (1150/1223-1246/1309) -39.8 W06-W12 N08-N13
      9 1572 1216 (N10W34 L=269)   627    145        3 (1500/1192-1647/1252)   0.0 W28-W39 N08-N12
     10 1765 1113 (N06W50 L=285)   474     98        3 (1737/1073-1813/1147) -23.2 W47-W54 N03-N08
     11  437 1040 (S05E38 L=197)   824    166        2 ( 415/ 960- 466/1123) -90.0 E39-E35 S10-N00
     12  835 1058 (S03E12 L=223)   241     58        1 ( 799/1042- 868/1078) -14.0 E13-E09 S04-S02
     13 1855  982 (S01W60 L=295)   432     58        1 (1836/ 948-1877/1010)  45.0 W57-W62 S02-N00
     14  148  882 (S15E66 L=169)  1324    284        5 ( 119/ 806- 177/ 971)  56.3 E70-E61 S18-S09
     15 1120  756 (S20W07 L=242)  1095    258        1 (1023/ 659-1254/ 837)  40.6 W01-W15 S26-S14
     16 1777  549 (S27W62 L=297)  1868    281        4 (1694/ 502-1843/ 622)  31.6 W53-W67 S30-S22
     17  941  548 (S34E03 L=232)   202     37        0 ( 926/ 528- 965/ 570) -18.4 E04-E01 S35-S32
     18  862  498 (S38E09 L=226)   288     95        0 ( 816/ 470- 908/ 530)  35.5 E12-E05 S40-S35
  

    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