Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  520 1377 (N11E39 L=352)   349    113        0 ( 506/1327- 538/1425) -45.0 E41-E36 N08-N13
      2  904 1299 (N13E12 L= 19)   348    155        1 ( 876/1226- 943/1377)  90.0 E13-E11 N08-N18
      3 1545  696 (S11W40 L= 72)   593    121        0 (1520/ 635-1564/ 755) -74.1 W37-W41 S14-S07
      4 1019  707 (S20W06 L= 37)   259    125        2 ( 960/ 687-1085/ 738)   0.0 W01-W09 S21-S19
      5 1118  617 (S24W14 L= 46)   343    123        0 (1059/ 589-1183/ 650)  35.5 W10-W17 S26-S20
      6 1260  547 (S25W25 L= 57)   677    240        1 (1154/ 489-1350/ 580)  -7.1 W17-W33 S27-S23
      7 1448  519 (S24W39 L= 71)   460     59        2 (1415/ 486-1476/ 548)  53.1 W37-W40 S25-S21
      8 1552  393 (S29W54 L= 86)  1859    368        1 (1413/ 333-1714/ 472)  36.9 W45-W66 S36-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