Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1702 1429 (N22W57 L=189)   351     46        1 (1677/1405-1723/1449) -45.0 W54-W59 N20-N23
      2 1805 1383 (N18W68 L=200)   836     74        3 (1767/1359-1838/1403)  -4.8 W62-W74 N15-N19
      3 1743 1369 (N18W59 L=191)   360     55        3 (1721/1337-1758/1397)  45.0 W56-W61 N15-N19
      4 1871 1112 (N00W66 L=198)   597     72        1 (1844/1077-1887/1150)  45.0 W62-W69 S01-N02
      5 1177 1071 (N07W10 L=142)   474    104        2 (1154/1031-1213/1129) -90.0 W08-W11 N04-N10
      6 1371  804 (S12W20 L=152)   377    160        0 (1346/ 747-1404/ 855)  63.4 W18-W22 S15-S08
      7 1767  741 (S21W53 L=185)  1992    242        6 (1679/ 664-1840/ 842)  31.0 W44-W60 S25-S15
      8 1180  745 (S13W07 L=139)   368    152        0 (1106/ 725-1231/ 769)  -7.1 W02-W10 S14-S11
      9 1575  518 (S34W38 L=170)   641     90        2 (1537/ 483-1615/ 554)  59.0 W34-W41 S36-S31
     10  397  506 (S24E54 L= 78)   721    105        3 ( 359/ 466- 447/ 536)   8.7 E61-E48 S27-S22
     11 1390  312 (S49W26 L=158)   479     52        2 (1365/ 287-1415/ 344) -50.2 W22-W28 S52-S45
  

    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