Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  939 1592 (N41E07 L=353)   283     41        0 ( 910/1576- 964/1608)   0.0 E09-E04 N39-N42
      2 1508 1379 (N25W35 L= 34)   424    107        1 (1488/1334-1536/1429)  80.5 W33-W36 N22-N28
      3  599 1274 (N19E29 L=331)   969    307        0 ( 528/1166- 685/1391) -43.0 E36-E21 N12-N26
      4 1222 1166 (N12W13 L= 12)   924    239        1 (1132/1094-1309/1245) -45.0 W07-W18 N07-N17
      5  999 1108 (N09E02 L=358)   307    132        0 ( 977/1050-1017/1166) -81.9 E02-W00 N05-N12
      6 1705 1045 (N03W47 L= 46)   409    103        0 (1683/1000-1721/1093)  80.5 W45-W48 N00-N06
      7 1091  795 (S11W04 L=  3)   680    209        3 (1031/ 750-1175/ 857)  29.1 W00-W09 S13-S06
      8  979  752 (S13E03 L=356)   224     56        0 ( 963/ 715- 996/ 783) -76.0 E03-E01 S15-S11
      9 1638  607 (S24W46 L= 45)   425    121        1 (1588/ 578-1687/ 651)  33.7 W41-W50 S26-S21
  

    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