Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1343 1589 (N35W05 L= 23)   372     60        0 (1316/1565-1373/1610)  59.0 W03-W06 N32-N37
      2  248 1478 (N04E70 L=309)  1101    143        2 ( 219/1446- 289/1502)  33.7 E75-E65 N01-N05
      3  129 1206 (S14E72 L=307)  2736    421        4 (  87/1092- 171/1326) -66.8 E82-E63 S20-S06
      4  351 1100 (S17E45 L=334)  1248    282        1 ( 288/ 973- 387/1240)  83.3 E48-E42 S25-S08
      5  634 1104 (S12E25 L=354)   556    139        3 ( 595/1043- 679/1150) -80.5 E27-E22 S14-S08
      6  780 1072 (S10E15 L=  4)   188     73        0 ( 739/1048- 818/1089)   0.0 E17-E12 S10-S09
      7 1911  988 (N17W63 L= 81)  1452    267        4 (1883/ 937-1957/1034)  81.9 W58-W69 N13-N20
      8 1715  759 (S02W51 L= 69)  1545    335        1 (1622/ 717-1813/ 846)   0.0 W40-W61 S04-N00
      9 1124  752 (S20W13 L= 32)   477    160        0 (1101/ 679-1149/ 834) -60.9 W09-W15 S24-S15
     10  781  707 (S31E07 L= 12)  1008    329        1 ( 610/ 667- 931/ 739)  15.3 E18-W04 S35-S28
     11 1524  426 (S26W52 L= 71)  2662    502        6 (1364/ 356-1715/ 519)  39.3 W43-W65 S34-S15
     12 1632  383 (S24W64 L= 83)  1677    172        3 (1584/ 333-1665/ 448) -63.4 W61-W67 S27-S19
  

    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