Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  407 1624 (N45E55 L=311)   497    151        0 ( 371/1603- 437/1634)   0.0 E58-E50 N43-N45
      2  534 1583 (N42E40 L=326)   425     77        2 ( 511/1561- 558/1602)  56.3 E42-E37 N40-N43
      3  330 1481 (N34E56 L=310)  1624    334        4 ( 289/1400- 379/1552) -37.6 E62-E48 N28-N38
      4 1468 1507 (N33W36 L= 41)   357     91        0 (1424/1483-1513/1532) -26.6 W32-W39 N31-N34
      5  898 1455 (N32E08 L=358)   376     94        0 ( 853/1438- 949/1470) -15.9 E10-E03 N30-N33
      6 1278 1405 (N27W19 L= 24)   520    169        0 (1209/1375-1369/1424) -15.3 W14-W25 N24-N28
      7  683 1278 (N21E22 L=343)  2736    744        9 ( 432/1063- 846/1423) -35.4 E42-E10 N07-N30
      8 1143 1381 (N26W09 L= 15)   299     54        1 (1106/1354-1177/1405)  31.0 W06-W11 N24-N27
      9  252 1187 (N15E57 L=309)  1734    342        6 ( 213/1069- 294/1315) -61.9 E61-E51 N07-N22
     10 1533 1201 (N13W34 L= 40)   643    183        0 (1457/1129-1606/1264) -45.0 W29-W39 N07-N17
     11  361 1055 (N07E45 L=320)   496    165        1 ( 311/ 999- 400/1126)  48.8 E49-E42 N03-N11
     12  632  616 (S21E28 L=337)  1624    559        3 ( 451/ 492- 797/ 764)  31.0 E45-E15 S29-S11
     13 1386  743 (S15W23 L= 28)   249     57        1 (1353/ 727-1415/ 762)  26.6 W20-W24 S15-S13
     14 1749  629 (S25W54 L= 60)   485     79        0 (1724/ 615-1777/ 645) -26.6 W51-W57 S26-S23
     15  947  531 (S27E07 L=358)   622    126        1 ( 889/ 504-1027/ 565)  11.3 E11-E01 S29-S25
     16  837  453 (S33E16 L=350)   437    117        1 ( 812/ 398- 864/ 490)  54.5 E18-E13 S37-S30
  

    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