Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1393 1713 (N43W49 L= 13)   425     77        0 (1363/1692-1423/1731)   0.0 W44-W52 N41-N44
      2  626 1434 (N37E24 L=301)   338     75        1 ( 588/1403- 663/1464) -33.7 E26-E20 N34-N39
      3 1551 1298 (N14W39 L=  3)   767    173        1 (1521/1233-1593/1375) -84.3 W36-W41 N09-N19
      4 1158 1228 (N16W11 L=335)   323    159        0 (1074/1199-1216/1263) -32.0 W06-W14 N13-N19
      5  413 1164 (N22E40 L=284)   382    123        1 ( 377/1124- 461/1203) -40.6 E43-E36 N18-N24
      6  828 1152 (N16E10 L=314)   259     52        0 ( 807/1123- 844/1179) -90.0 E11-E09 N14-N18
      7 1691  701 (S25W43 L=  7)   763    182        1 (1641/ 630-1739/ 776)  45.0 W38-W47 S29-S21
      8  496  406 (S28E52 L=272)  2217    275        5 ( 463/ 291- 538/ 528)  32.3 E72-E42 S37-S18
      9  965  324 (S40E17 L=307)   308     64        0 ( 934/ 308- 999/ 340)   0.0 E20-E14 S41-S38
  

    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