Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  137 1222 (S12E74 L=291)  1042    153        0 ( 113/1170- 167/1273)  66.8 E76-E70 S15-S08
      2 1672 1143 (N20W38 L= 43)   615    127        1 (1629/1123-1725/1169)   7.1 W34-W42 N18-N21
      3  320 1058 (S18E47 L=318)   395    126        0 ( 272/1024- 363/1077)  35.5 E50-E43 S21-S16
      4 1518  983 (N07W30 L= 35)   286     83        0 (1485/ 934-1542/1023)  80.5 W29-W31 N03-N09
      5  879  764 (S22E02 L=  3)   587    161        1 ( 802/ 727- 939/ 823)  60.9 E05-W00 S26-S17
      6 1234  710 (S16W21 L= 26)   332    114        0 (1168/ 688-1288/ 735)  35.5 W16-W23 S18-S13
      7  604  654 (S37E18 L=347)   230     46        1 ( 587/ 627- 631/ 678) -14.0 E19-E15 S38-S35
      8 1214  565 (S25W25 L= 30)   311    112        0 (1153/ 552-1270/ 579)  20.6 W20-W28 S26-S22
      9  835  500 (S40W03 L=  8)   898    175        3 ( 784/ 427- 903/ 571) -30.3 E02-W10 S43-S36
     10 1020  436 (S39W19 L= 24)   433    121        0 ( 952/ 415-1073/ 459)  18.4 W13-W22 S40-S37
     11 1539  419 (S24W54 L= 59)   639    172        1 (1481/ 401-1607/ 431)  26.6 W50-W60 S26-S21
  

    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