Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1583 (N27E31 L=228)  2237    366        7 ( 669/1429- 826/1720) -35.3 E45-E21 N17-N35
      2 1195 1705 (N47E05 L=255)   498    105        2 (1137/1691-1247/1720)  29.1 E09-W00 N44-N49
      3  585 1385 (N11E35 L=224)   261     77        1 ( 560/1369- 623/1406)   0.0 E37-E32 N09-N12
      4 1681 1289 (N28W41 L=300)   553    112        2 (1634/1254-1734/1318)   0.0 W36-W46 N27-N29
      5 1266 1213 (N14W11 L=270)   507    149        4 (1233/1147-1297/1258) -74.1 W08-W12 N10-N17
      6 1268 1055 (N05W14 L=273)   200     81        0 (1229/1035-1305/1068)  36.9 W11-W15 N03-N06
      7 1052  853 (S11W05 L=265)   532    103        3 ( 996/ 832-1104/ 888)   0.0 W01-W08 S12-S08
      8 1108  703 (S19W12 L=271)   444    197        0 (1022/ 679-1192/ 732)  33.7 W07-W16 S21-S15
      9  474  387 (S58E40 L=220)  1157    156        2 ( 405/ 364- 545/ 411)   9.5 E54-E27 S61-S54
     10 1015  371 (S43W17 L=277)   494    142        0 ( 949/ 350-1076/ 407)  40.6 W12-W20 S46-S38
     11 1334  350 (S36W43 L=302)   863    155        2 (1264/ 327-1393/ 379)  21.8 W35-W47 S38-S33
  

    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