Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1284 (N02E45 L= 52)   284     41        0 ( 385/1264- 434/1302)  26.6 E47-E43 N00-N03
      2 1190 1218 (N13W06 L=103)   511    163        0 (1137/1170-1257/1290) -29.1 W01-W10 N11-N16
      3  576  984 (S13E27 L= 70)   325    121        0 ( 561/ 929- 590/1051) -74.1 E28-E25 S15-S08
      4 1171  933 (S04W10 L=107)   511    196        0 (1104/ 841-1206/1024) -41.6 W04-W13 S08-N00
      5 1753  944 (N08W50 L=147)   532    156        1 (1700/ 912-1800/ 972)  -5.7 W45-W55 N07-N09
      6  846  656 (S27E04 L= 93)  1357    563        0 ( 638/ 533-1104/ 755)  42.7 E17-W10 S40-S16
      7 1365  623 (S18W30 L=127)   290    102        0 (1321/ 579-1414/ 648)  45.0 W28-W33 S21-S16
      8 1467  613 (S17W38 L=134)   300     44        0 (1446/ 588-1489/ 637) -76.0 W35-W38 S18-S14
      9 1526  561 (S19W44 L=141)   309     72        1 (1499/ 539-1549/ 593)  76.0 W42-W45 S20-S16
     10  885  509 (S36W02 L= 99)   277     62        1 ( 860/ 483- 914/ 537)  68.2 W00-W03 S38-S33
     11 1549  464 (S24W51 L=147)   419    192        0 (1490/ 453-1620/ 480)  14.0 W45-W57 S25-S22
     12 1444  278 (S39W55 L=152)  1100    194        1 (1376/ 259-1512/ 304)  29.1 W49-W59 S40-S35
  

    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