Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1143 1599 (N33W14 L=271)   522    100        2 (1100/1565-1180/1627) -11.3 W10-W16 N30-N35
      2 1407 1203 (N06W25 L=281)   755    189        3 (1323/1155-1495/1238) -20.0 W19-W30 N02-N08
      3 1796  848 (S18W53 L=309)  1013    172        3 (1755/ 776-1836/ 914)  76.0 W49-W56 S22-S13
      4  290  858 (S04E51 L=205)   747    116        2 ( 261/ 809- 329/ 909) -71.6 E53-E47 S07-S01
      5 1623  697 (S26W39 L=295)   328     87        2 (1582/ 670-1665/ 724)   8.1 W35-W42 S27-S24
      6 1738  660 (S29W51 L=307)   892    207        2 (1649/ 619-1828/ 687)   6.0 W42-W62 S31-S27
      7  385  598 (S21E50 L=206)   587    161        2 ( 357/ 523- 405/ 669)  60.9 E54-E47 S25-S16
      8  974  507 (S33E09 L=248)  1077    392        3 ( 800/ 433-1135/ 551)   4.8 E22-W02 S36-S30
      9 1408  462 (S40W24 L=280)   915    270        1 (1296/ 434-1529/ 530)  14.7 W14-W33 S41-S36
     10  381  389 (S34E65 L=191)   622     79        0 ( 360/ 355- 407/ 420)  14.9 E74-E59 S36-S32
     11  661  389 (S38E36 L=220)   777    185        4 ( 575/ 359- 747/ 416)  -3.8 E44-E28 S38-S36
     12 1523  339 (S52W41 L=297)  1413    181        3 (1450/ 312-1610/ 366)   0.0 W31-W54 S53-S48
  

    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