Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  549 1661 (N58E36 L=  9)  2926    296        5 ( 427/1590- 682/1723)   6.5 E55-E18 N52-N61
      2  854 1571 (N44E02 L= 44)   661    146        1 ( 785/1549- 926/1591) -26.6 E07-W03 N41-N46
      3 1305 1544 (N33W31 L= 77)   258     47        0 (1279/1521-1323/1565) -90.0 W29-W32 N31-N34
      4 1073 1533 (N37W14 L= 59)   430    150        1 (1003/1515-1142/1550)  -5.2 W08-W19 N35-N38
      5 1189 1318 (N21W16 L= 61)   784    202        1 (1119/1255-1269/1345)  -5.7 W10-W20 N18-N23
      6  263 1270 (N31E57 L=348)   471     90        0 ( 225/1253- 316/1286)  -4.8 E62-E50 N30-N32
      7  776 1102 (N15E14 L= 32)   644    174        1 ( 735/1022- 801/1179) -78.7 E15-E12 N10-N20
      8  289 1059 (N18E52 L=353)   909    186        0 ( 265/ 966- 316/1124) -73.3 E54-E49 N12-N22
      9 1822  863 (S19W56 L=101)  1237    221        0 (1773/ 777-1869/ 971)  39.3 W50-W62 S23-S13
     10 1583  714 (S23W32 L= 78)  1212    249        2 (1486/ 654-1642/ 811)  28.3 W24-W37 S25-S17
     11 1137  680 (S16W01 L= 47)   627    152        2 (1108/ 611-1181/ 743)  63.4 E01-W04 S20-S12
     12 1649  653 (S29W38 L= 84)   297     86        0 (1620/ 630-1690/ 683)  15.9 W35-W42 S29-S27
     13 1221  525 (S28W04 L= 50)   576    113        2 (1180/ 480-1260/ 572) -63.4 W02-W06 S32-S24
  

    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