Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  690 1567 (N31E32 L=287)   723    136        1 ( 665/1512- 721/1625) -63.4 E34-E28 N27-N35
      2  540 1405 (N19E38 L=281)   341    122        0 ( 528/1353- 558/1464) -74.1 E39-E35 N15-N22
      3 1720 1328 (N28W49 L=  9)   738    194        0 (1666/1242-1773/1397) -48.4 W45-W54 N22-N32
      4 1274 1166 (N12W14 L=333)   714    294        1 (1209/1038-1315/1290) -59.0 W08-W17 N04-N19
      5 1562 1101 (N11W34 L=353)   761    177        0 (1536/1013-1580/1188)  84.8 W32-W35 N05-N16
      6  393  940 (S12E41 L=278)   408     96        1 ( 370/ 897- 417/ 983)  80.5 E43-E39 S14-S08
      7  198  876 (S18E62 L=257)   986      0       -1 ( 151/ 843- 232/ 908) -10.3 E69-E57 S19-S16
      8 1560  495 (S26W48 L=  7)  2455    383        1 (1413/ 420-1720/ 576)  29.5 W37-W61 S33-S19
      9  893  511 (S33E03 L=317)   352    106        0 ( 882/ 455- 903/ 567) -76.0 E03-E01 S37-S29
  

    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