Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1537 1610 (N32W51 L= 33)   414    150        0 (1495/1582-1563/1644)  12.5 W46-W55 N30-N33
      2 1259 1489 (N31W25 L=  7)   723    213        1 (1176/1440-1361/1529)  -8.1 W19-W33 N28-N32
      3 1748 1421 (N16W61 L= 43)  1336    211        1 (1692/1394-1811/1451) -28.6 W55-W67 N12-N18
      4 1335  951 (S03W17 L=360)  1923    467        5 (1274/ 741-1378/1160)  73.7 W10-W20 S14-N10
      5  718 1095 (N15E18 L=325)   505    160        1 ( 649/1056- 787/1145) -36.9 E21-E13 N11-N19
      6 1612  955 (S08W36 L= 19)   697    205        3 (1549/ 885-1690/1008) -48.0 W32-W41 S14-S03
      7 1515  671 (S24W27 L=  9)   884    182        0 (1482/ 599-1536/ 770)  65.6 W23-W28 S28-S17
      8 1120  385 (S37E05 L=337)  1439    307        2 (1033/ 304-1188/ 503)  40.6 E13-W01 S42-S28
  

    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