Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1433 1751 (N36W55 L=137)  1312    126        3 (1377/1712-1479/1780) -50.2 W49-W61 N32-N38
      2 1720 1517 (N15W63 L=145)   427     47        0 (1693/1498-1748/1532) -31.0 W59-W65 N13-N16
      3  757 1454 (N34E08 L= 74)  1660    348        0 ( 615/1371- 894/1513) -33.7 E18-W00 N25-N39
      4 1390 1453 (N19W33 L=115)   593    137        0 (1344/1403-1444/1508)  16.7 W28-W38 N16-N21
      5  242 1337 (N37E56 L= 26)  2019    172        3 ( 212/1260- 280/1424) -42.7 E64-E49 N31-N43
      6 1229 1385 (N19W20 L=102)   538    136        0 (1156/1362-1294/1409) -12.5 W15-W24 N17-N20
      7  227 1016 (N17E56 L= 26)  1535    267        2 ( 196/ 929- 284/1147) -77.9 E59-E50 N11-N25
      8 1888 1040 (S16W63 L=145)   625     85        0 (1860/1007-1911/1083)  29.7 W58-W67 S17-S13
      9 1065  785 (S12E03 L= 79)   585    194        1 (1034/ 698-1089/ 886)  63.4 E05-W00 S16-S06
     10 1572  748 (S26W29 L=111)   321     67        0 (1558/ 706-1594/ 778) -90.0 W28-W31 S28-S23
     11  798  602 (S17E23 L= 59)   304     59        1 ( 774/ 578- 831/ 631) -63.4 E24-E21 S19-S15
     12 1488  576 (S35W21 L=103)   257     48        0 (1465/ 555-1515/ 592)  14.0 W19-W23 S36-S33
     13 1380  510 (S37W12 L= 94)   239     49        0 (1353/ 495-1404/ 531)  11.3 W09-W14 S37-S35
     14  336  468 (S17E66 L= 16)  1374    139        5 ( 295/ 409- 367/ 528)  32.5 E73-E60 S20-S13
     15  477  417 (S22E55 L= 27)  1129    191        2 ( 424/ 336- 518/ 500)  63.4 E60-E51 S28-S16
     16  649  336 (S31E46 L= 36)   395     91        0 ( 608/ 309- 695/ 361)   5.7 E51-E41 S32-S29
  

    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