Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1324 1534 (N35W06 L= 77)  1307    278        4 (1241/1413-1372/1640) -47.3 E02-W11 N27-N41
      2  406 1488 (N10E54 L= 16)   801    126        3 ( 378/1427- 431/1546) -58.0 E57-E50 N05-N13
      3  707 1444 (N14E30 L= 41)   233    107        0 ( 690/1400- 740/1483)  90.0 E30-E28 N10-N16
      4  225 1414 (N02E71 L=359)   702     57        1 ( 200/1390- 249/1442)  26.6 E75-E68 N00-N03
      5 1629 1155 (N19W35 L=105)  1176    233        1 (1575/1058-1669/1265)  85.9 W31-W37 N13-N27
      6 1848 1031 (N19W58 L=129)  1612    179        2 (1808/ 954-1877/1121) -56.3 W53-W63 N14-N23
      7 1330 1012 (N03W18 L= 88)   423     99        0 (1280/ 985-1368/1045)  45.0 W15-W19 N01-N06
      8  199  999 (S23E61 L= 10)   410     78        2 ( 179/ 956- 221/1030)  53.1 E63-E57 S25-S21
      9  552  697 (S34E23 L= 47)   957    197        0 ( 474/ 668- 655/ 741)   7.6 E30-E15 S35-S31
     10  914  702 (S25W01 L= 72)   284     65        0 ( 878/ 684- 951/ 720)  33.7 W00-W03 S26-S22
     11 1700  580 (S11W60 L=130)   768    168        3 (1652/ 521-1751/ 659)  78.7 W56-W61 S15-S05
     12 1595  474 (S20W57 L=127)  1428    342        4 (1521/ 366-1666/ 596) -76.0 W51-W66 S27-S11
  

    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