Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1581 1654 (N46W54 L=255)  1535    253        3 (1489/1614-1689/1694)  -8.4 W42-W71 N42-N49
      2 1038 1558 (N37E01 L=201)   369     56        0 (1006/1538-1079/1577)  11.3 E03-W02 N35-N38
      3  334 1497 (N29E60 L=141)   767    116        0 ( 308/1452- 379/1524)  63.4 E65-E55 N26-N31
      4 1103 1483 (N32W04 L=206)   358     86        0 (1076/1436-1127/1515) -59.0 W02-W05 N29-N34
      5 1100 1344 (N23W04 L=206)   382    110        0 (1076/1301-1149/1375)  45.0 W02-W07 N20-N24
      6  186 1172 (N08E65 L=137)   709    127        0 ( 172/1111- 205/1235) -49.4 E67-E61 N04-N11
      7  774 1191 (N12E16 L=185)   235     83        0 ( 737/1165- 813/1215) -31.0 E18-E13 N10-N13
      8 1499 1099 (N08W30 L=232)  1233    291        2 (1416/ 965-1558/1196) -57.3 W24-W34 N00-N14
      9 1654 1063 (N06W42 L=244)   415     80        1 (1626/1037-1687/1107) -45.0 W39-W44 N04-N08
     10 1733  816 (S09W50 L=252)  1184    173        2 (1690/ 743-1769/ 899)  68.2 W47-W53 S13-S03
     11  203  842 (S12E62 L=139)   543     68        1 ( 178/ 821- 229/ 863)   0.0 E65-E58 S13-S10
     12 1117  764 (S13W06 L=208)   198     57        0 (1095/ 739-1144/ 794) -45.0 W05-W08 S14-S10
     13  193  733 (S19E67 L=135)   800     83        0 ( 170/ 693- 214/ 774)  33.7 E74-E62 S22-S16
     14  702  603 (S25E21 L=181)  1067    225        2 ( 626/ 567- 792/ 657)  -5.2 E25-E14 S27-S21
     15 1642  595 (S23W47 L=249)  1503    178        2 (1582/ 523-1698/ 657)  56.3 W42-W51 S28-S19
     16  499  594 (S26E37 L=165)   365     98        0 ( 453/ 578- 534/ 623)  15.9 E40-E33 S27-S24
     17 1192  220 (S55W22 L=224)  2271    254        5 (1102/ 189-1314/ 250)   7.4 W12-W36 S59-S51
  

    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