Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  432 1702 (N32E57 L=202)   668    141        1 ( 417/1674- 452/1736) -35.5 E60-E53 N29-N35
      2 1163 1473 (N23W03 L=261)   746    178        1 (1100/1435-1219/1523)  26.6 W00-W06 N20-N26
      3  204 1360 (N08E65 L=193)  2330    326        3 ( 126/1303- 268/1400)  15.5 E77-E57 N04-N10
      4 1395 1228 (N11W20 L=278)  1258    278        4 (1339/1139-1446/1307) -48.8 W15-W23 N06-N15
      5  966 1230 (N06E06 L=252)   214     60        0 ( 925/1212- 996/1254)  26.6 E08-E04 N04-N07
      6  255 1180 (S03E55 L=204)   576    149        1 ( 224/1126- 300/1240) -63.4 E58-E50 S06-N00
      7  166 1088 (S09E63 L=195)   727    133        2 ( 138/1042- 205/1127)  39.8 E66-E58 S11-S06
      8  703  839 (S20E18 L=240)   206     68        1 ( 681/ 814- 742/ 868) -45.0 E19-E15 S21-S18
      9 1494  803 (S12W32 L=290)   455    124        1 (1429/ 778-1547/ 830)  23.2 W28-W35 S14-S10
     10 1354  755 (S17W24 L=282)   242     47        0 (1320/ 734-1381/ 772) -11.3 W21-W26 S18-S16
     11 1510  403 (S36W47 L=306)   317    118        0 (1466/ 389-1562/ 413)  23.2 W44-W51 S37-S34
     12  866  203 (S66W01 L=259)   631     89        1 ( 838/ 169- 893/ 227) -36.9 E03-W05 S69-S63
  

    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